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Abstract During the past two decades, it has been demonstrated that chiral effec- 
tive field theory represents a powerful tool to deal with nuclear forces in a systematic 
and model-independent way. Two-, three-, and four-nucleon forces have been derived 
up to next-to-next-to-next-to-leading order (N'^LO) and (partially) applied in nuclear 
few- and many-body systems — with, in general, a good deal of success. This may sug- 
gest that we are finally done with the nuclear force problem; but that would be too 
optimistic. There are still some pretty basic open issues that have been swept under 
rug and, finally, need our full attention, like the proper renormalization of the two- 
nucleon potential. Moreover, the order-by-order convergence of the many-body force 
contributions is at best obscure at this time. 

Keywords Effective field theory ■ Chiral perturbation theory • Renormalization • 
Nuclear two- and many-body forces 



1 Introduction 



The theory of nuclear forces has a long history (cf. Table [TJ. Based upon the seminal 
idea by Yukawa [1], first field-theoretic attempts to derive the nucleon-nucleon {NN) 
interaction focused on pion exchange. While the one-pion exchange turned out to be 
very useful in explaining A'^A'^ scattering data and, in particular, the properties of the 
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Table 1 The Theory of Nuclear Forces: Eight Decades of Struggle 



1935 


Yukawa: Meson Theory 




The "Pion Theories". 


1950's 


One-pion exchange: good; 




Multi-pion exchange: disaster. 


1960's 


Many pions = multi-pion resonances: cr, p, lj, ... 




The One-Boson-Exchange Model: success. 


1970's 


Refinement of meson theory: Sophisticated 27r exchange models; 




Partovi-Lomon, Stony Brook, Paris, Bonn. 


1980's 


Nuclear physicists discover QCD: 




Quark Models. 


1990's 


Nuclear physicists discover EFT; Weinberg, van Kolck, . . . 


and beyond 


Back to Yukawa's Meson (Pion) Theory! 




But, constrained by Chiral Symmetry. 



deuteron [2], multi-pion exchange was beset with serious ambiguities !3! that could not 
be resolved in a satisfactory way. Thus, the "pion theories" of the 1950s are generally 
judged as failures — for reasons we understand today: pion dynamics is constrained by 
chiral symmetry, a crucial point that was unknown in the 1950s. 

Historically, the experimental discovery of heavy mesons in the early 1960s saved the 
situation. The one-boson-exchange (OBE) model emerged, which still today is the most 
economical and quantitative phenomenology for describing the A'^A'^ interaction [1]. 
The weak point of this model, however, is the scalar-isoscalar "sigma" or "epsilon" 
boson, for which empirical evidence remains controversial. Since this boson is associated 
with the correlated (or resonant) exchange of two pions, a vast theoretical effort was 
launched to derive the 27r-exchange contribution of the nuclear force, which creates the 
intermediate range attraction. During this effort, which occupied more than a decade, 
dispersion theory (Stony Brook and Paris potentials) as well as field theory (Partovi- 
Lomon model, Bonn potential l4]|5]) were invoked. 

The nuclear force problem appeared to be solved; however, with the discovery of 
quantum chromo-dynamics (QCD), all "meson theories" were relegated to the status 
of models and the attempts to derive the nuclear force had to start all over again. 

The problem with a derivation of nuclear forces from QCD is that this theory 
is non-perturbative in the low-energy regime characteristic of nuclear physics, which 
makes direct solutions very difficult. Therefore, during the first round of new attempts, 
QCD-inspired quark models became popular. The positive aspect of these models is 
that they try to explain hadron structure and hadron-hadron interactions on an equal 
footing and, indeed, some of the gross features of the A'^A'^ interaction are explained 
successfully. However, on a critical note, it must be pointed out that these quark- 
based approaches are nothing but another set of models and, thus, do not represent 
fundamental progress. For the purpose of describing hadron-hadron interactions, one 
may equally well stay with the simpler and much more quantitative meson models. 

A major breakthrough occurred when the concept of an effective field theory (EFT) 
was introduced and applied to low-energy QCD. As outlined by Weinberg in a seminal 
paper j6], one has to write down the most general Lagrangian consistent with the 
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assumed symmetry principles, particularly the (broken) chiral symmetry of QCD. At 
low energy, the effective degrees of freedom are pious (the Goldstone bosons of the 
broken symmetry) and nucleons rather than quarks and gluous; heavy mesons and 
nucleon resonances are "integrated out" . So, the circle of history is closing and we are 
back to Yukawa's meson (pion) theory, except that we have finally learned how to deal 
with it: broken chiral symmetry is a crucial constraint that generates and controls the 
dynamics and establishes a clear connection with the underlying theory, QCD. 

The past 15 years have seen great progress in applying chiral perturbation theory 
(ChPT) to nuclear forces [7ll8ll9ll 1011 11 IIT21IT3UT4] . As a resuh, nucleon- nucleon (NN) 
potentials of high precision have been constructed, which are based upon ChPT car- 
ried to next-to-next-to-next-to- leading order (N'^LO) 12,13,14 , and applied in nuclear 
structure calculations with great success. 

However, in spite of this progress, we are not done. Due to the complexity of the 
nuclear force issue, there are still many subtle and not so subtle open problems. We 
will not list and discuss all of them, but instead just focus on the two open issues, 
which we perceive as the most important ones: 

— The proper renormalization of chiral nuclear potentials and 

— Subleading chiral few-nucleon forces. 



2 Renormalization of chiral nuclear forces 

2.1 The chiral NN potential 

In terms of naive dimensional analysis or "Weinberg counting" , the various orders of 
the irreducible graphs which define the chiral TV TV potential are given by: 

= + vi:^ (1) 

Vnlo = VLo + V^P + Vl^^ + V^l^ (2) 

^NNLO = ^NLO + V}^^ + V^^^ (3) 

V^N3LO = VkNLO + V;^^ + + V^t^ + V^f (4) 

where the superscript denotes the order u of the low-momentum expansion. LO stands 
for leading order, NLO for next-to- leading order, etc.. Contact potentials carry the 
subscript "ct" and pion-exchange potentials can be identified by an obvious subscript. 
For more details concerning the above potentials, see ref. |14| . 



2.2 Nonperturbative renormalization of the TV TV potential 

The two-nucleon system is characterized by large scattering lengths and shallow (quasi) 
bound states which require a nonperturbative treatment. Following Weinberg's pre- 
scription |7], this is accomplished by inserting the potential V into the Lippmann- 
Schwinger (LS) equation: 

T(P P) = ^^(p ', P) + / dV V(p ', p ") ^% . T(p ", p) , (5) 

J — p" + le 



where Mj^ denotes the nucleon mass. 
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In general, the integral in the LS equation is divergent and needs to be regularized. 
One way to achieve this is by multiplying V with a regulator function 

V{p ',p) ^ V{p ',p) e-('''/^)- ^-(P/Af- . (6) 

Typical choices for the cutoff parameter A that appears in the regulator are A « 
0.5 GeV <Ax^l GeV. 

In field theories, divergent integrals are not uncommon and methods have been 
designed to deal with them. One regulates the integrals and then removes the depen- 
dence on the regularization parameters (scales, cutoffs) by "renormalization" . In the 
end, the theory and its predictions do not depend on cutoffs or renormalization scales. 
So-called renormalizable quantum field theories, like QED, have essentially one set of 
prescriptions that takes care of renormalization through all orders. In contrast, EFTs 
are renormalized by "counter terms" (contact terms) that are introduced order by order 
in increasing numbers. 

Naively, the most perfect renormalization procedure is the one where the cutoff 
parameter A is taken to infinity while stable and quantitative results are maintained 
through the adjustment of counter terms. This was accomplished at LO in the work 
by Nogga et al ^5j. At NNLO, the infinite-cutoff renormalization procedure has been 
investigated in [TS] for partial waves with total angular momentum J < 1 and in [T7] for 
all partial waves with J < 5. However, for a quantitative chiral A'' A'' potential one needs 
to advance all the way to N^LO. At N^LO, the state was considered in Ref. |18| . 
and all states up to J = 6 were investigated in Ref. [19]. From all of these works, it is 
evident that no counter term is effective in partial- waves with short-range repulsion and 
only a single counter term can constructively be used in partial-waves with short-range 
attraction. Thus, for the yl — ^ oo renormalization prescription, even at N'^LO, there 
exists either one or no counter term per partial-wave state. This is inconsistent with 
any reasonable power-counting scheme and prevents an order-by-order improvement of 
the predictions. 

To summarize: In the infinite-cutoff renormalization scheme, the potential is ad- 
mitted up to unlimited momenta. However, the EFT this potential is derived from 
has validity only for momenta smaller than the chiral symmetry breaking scale yl^ « 
GeV. The lack of order-by-order convergence and discrepancies in lower partial-waves 
demonstrate that the potential should not be used beyond the limits of the effective 
theory [19J (see Ref. |20| for a related discussion). The conclusion then is that cutoffs 
should be limited to A < A^ (but see also Ref. [21] and J. Gegelia's contribution to 
this conference). 

Crucial for an EFT are regulator independence (within the range of validity of the 
EFT) and a power counting scheme that allows for order-by-order improvement with 
decreasing truncation error. The purpose of renormalization is to achieve this regulator 
independence while maintaining a functional power counting scheme. 

Thus, in the spirit of Lepage [22], the cutoff independence should be examined for 
cutoffs below the hard scale and not beyond. Ranges of cutoff independence within the 
theoretical error are to be identified using 'Lepage plots' 22 . Recently, we have started 
a systematic investigation of this kind. In our work, we quantify the error of the pre- 
dictions by calculating the /datum for the reproduction of the neutron-proton (np) 
elastic scattering data as a function of the cutoff parameter A of the regulator function 
Eq. ([6|. We have investigated the predictions by chiral np potentials at order NLO 
and NNLO applying Weinberg counting for the counter terms (NN contact terms). 
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Fig. 1 x^/datum for the reproduction of the np data in the energy range 35-125 MeV (upper 
frame) and 125-183 MeV (lower frame) as a function of the cutoff parameter A of the regulator 
function Eq. (jcjl. The (black) dashed curves show the /datum achieved with np potentials 
constructed at order NLO and the (red) solid curves are for NNLO. 



We show our results for the energy range 35-125 MeV in the upper frame of Fig. [T] 
and for 125-183 MeV in the lower frame. It is seen that the reproduction of the np 
data at these energies is generally poor at NLO, while at NNLO the /datum assumes 
acceptable values (a clear demonstration of order-by-order improvement). Moreover, at 
NNLO one observes "plateaus" of constant low ^or cutoff parameters ranging from 
450 to 850 MeV. This may be perceived as cutoff independence (and, thus, successful 
renormalization) in the relevant range of cutoff parameters. 



3 Few-nucleon forces and what is missing 

We will now discuss the other issue we perceive as unfinished and important, namely, 
subleading chiral few-nucleon forces. 

Nuclear three-body forces in ChPT were initially discussed by Weinberg [8]. The 
three-nucleon force (3NF) at NNLO, was derived by van Kolck [10] and applied, for 
the first time, in nucleon-deuteron scattering by Epelbaum et al. [23]. The leading 
4NF (at N^LO) was constructed by Epelbaum [22 and found to contribute in the 
order of 0.1 MeV to the ^He binding energy (total ^He binding energy: 28.3 MeV) 
in a preliminary calculation , confirming the traditional assumption that 4NF are 
essentially negligible. Therefore, the focus is on 3NF. 

For the order of a 3NF, we have 



i 



(7) 
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where L denotes the number of loops and Ai is the vertex index. We will use this 
equation to analyze 3NF contributions order by order. The first non-vanishing 3NF 
occurs at = 3 (NNLO), which is obtained when there are no loops (L = 0) and 

Ai = 1, i.e., Ai = 1 for one vertex while Ai — for all other vertices. There are 
three topologies which fulfill this condition, known as the two-pion exchange (2PE), 
one-pion exchange (fPE), and contact graphs. 

The 3NF at NNLO has been applied in calculations of few-nucleon reactions j26| . 
structure of light- and medium-mass nuclei ^27.,28, 29,.30, 31 ji32j , and nuclear and neu- 
tron matter [33II34||551|36| with a great deal of success. However, the famous 'Ay puz- 
zle' of nucleon-deuteron scattering [23l and the analogous problem with the analyzing 
power in p-^He scattering [37j is not resolved. Furthermore, the spectra of light nuclei 
leave room for improvement [27] . Since we are dealing with a perturbation theory, it is 
natural to turn to the next order when looking for improvements. 

The next order is N'^LO, where we have loop and tree diagrams. For the loops, we 
have L — 1 and, therefore, all Ai have to be zero to ensure = 4. Thus, these one-loop 
3NF diagrams can include only leading order vertices, the parameters of which are 
fixed from nN and A^A'' analysis. One sub-group of these diagrams (the 2PE graphs) 
has been calculated by Ishikawa and Robilotta |38) . and the other topologies have been 
evaluated by the Bochum-Bonn group ^39.40]. The N^LO 2PE 3NF has been applied in 
the calculation of nucleon-deuteron observables in Ref. [38] causing little impact. Very 
recently, the long-range part of the chiral N^LO 3NF has been tested in the triton [41] 
and in three- nucleon scattering ^42^ yielding only moderate effects. The long- and short- 
range parts of this force have been used in neutron matter calculations (together with 
the N^LO 4NF) producing relatively large contributions from the 3NF (4^. Thus, the 
ultimate assessment of the N'^LO 3NF is still outstanding and will require more few- 
and many-body applications. 

In the meantime, it is of interest to take already a look at the next order of 3NFs, 
which is N^LO or = 5 (of the Z\-less theory to which the present discussion is 
restricted because of lack of space). The loop contributions that occur at this order are 
obtained by replacing in the N'^LO loops one vertex by a Zlj = 1 vertex (with LEG Cj), 
which is why these loops may be more sizable than the N'^LO loops. The 2PE topology 
has already been evaluated [44] and turns out to be of modest size; moreover, it can 
be handled in a practical way by summing it up together with the 2PE topologies at 
NNLO and N^LO [44]. However, there are four more loop topologies, which are very 
involved and that have not been worked out yet. Finally, a tree topology at N'^LO 
provides a new set of 3N contact interactions, which have recently been derived by the 
Pisa group [15]. Contact terms are typically simple (as compared to loop diagrams) 
and their coefficients are unconstrained (except for naturalness). Therefore, it would 
be an attractive project to test some terms (tn particular, the spin-orbit terms) of the 
N^LO contact 3NF in calculations of few-body reactions (specifically, the p-d and 
p-^He Ay) and spectra of light nuclei. 

4 Conclusions and Outlook 

The past 15 years have seen great progress in our understanding of nuclear forces in 
terms of low-energy QCD. Key to this development was the realization that low-energy 
QCD is equivalent to an effective field theory (EFT) which allows for a perturbative 
expansion that has become known as chiral perturbation theory (ChPT). In this frame- 
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work, two- and many-body forces emerge on an equal footing and the empirical fact 
that nuclear many-body forces are substantially weaker then the two-nucleon force is 
explained automatically. 

In spite of the great progress and success of the past 15 years, there are still some 
unresolved issues. One problem is the proper renormalization of the chiral two- and 
many-nucleon potentials, where systematic investigations are already under way (of. 
Sec. 2). 

The other unfinished business are the few-nucleon forces beyond NNLO ("sub- 
leading few-nucleon forces") which arc needed to hopefully resolve some important 
outstanding nuclear structure problems. At orders N^LO and N^LO very many new 
3NF structures appear, some of which have already been tested. However, in view of 
the multitude of 3NF topologies it will take a while until we will have a proper overview 
of impact and convergence of these contributions. 

If the open issues discussed in this paper will be resolved within the next few years, 
then, after 80 years of desperate struggle, we may finally claim that the nuclear force 
problem is essentially under control. The greatest beneficiaries of such progress will be 
the fields of exact few-nucleon calculations and ab initio nuclear structure physics. 
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